Controlled Dicke Subradiance from a Large Cloud of Two-Level Systems 
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Dicke superradiance has been observed in many systems and is based on constructive interferences 
between many scattered waves. The counterpart of this enhanced dynamics, subradiance, is a 
destructive interference effect leading to the partial trapping of light in the system. In contrast to 
the robust superradiance, subradiant states are fragile and spurious decoherence phenomena hitherto 
obstructed the observation of such metastable states. We show that a dilute cloud of cold atoms is 
an ideal system to look for subradiance in free space and study various mechanisms to control this 
subradiance. 
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Interferences between scattered waves by many parti- 
cles can give rise to a large variety of phenomena, includ- 
ing collective effects such as Bragg or Mie scattering, well 
known in the context of classical optics [T] . More intrigu- 
ing situations can arise in the mesoscopic regime where 
interferences are at the origin of coherent backscattering 
and Anderson localization [2] . Interest in trapping waves 
in disordered media has spurred the efforts in multiple 
scattering of light, as photons seem to be an ideal can- 
didate for non-interacting waves. Several mechanisms 
allowing to increase the time an optical excitation can 
stay in a system are known to exist and are based on 
different physical phenomena [31 H] . Multiple scattering 
of light or radiation trapping for instance allows for in- 
creased photon trapping in the absence of interferences 
[5]. Anderson localization is another possibility where 
interferences form localized states with exponentially de- 
creasing coupling to the environment and a dense sam- 
ple fulfilling the Ioffe-Regel criterion [6] is assumed to 
be required. This letter addresses another fundamental 
mechanism leading to 'long-lived' modes of excitation: 
the study of subradiance in dilute clouds of cold atoms 
i.e. the trapping of light by destructive interferences in 
the single scattering regime. This mechanism is based 
on the pioneering work by Dicke who studied enhanced 
decay rates in small and large samples [TJ. Typically 
Dicke states are considered for an assembly of N two- 
level systems, realized e.g. by atoms [S] or quantum dots 
[9j. Different regimes can be studied, using either an 
initially fully inverted system with TV photons stored by 
the N atoms or an initial state in which the whole sys- 
tem shares a single excitation [TUl [TT]. These systems 
have attracted an increasing attention in the context of 
quantum information science |12H14j . where the accessi- 
ble Hilbert space can be restricted to single excitations 
using e.g. the Rydberg blockade [T5HT5] . 

In this letter, we will show how it is possible to under- 
stand and control the coupling of light into metastable 
subradiant states, illustrating that large dilute clouds of 
cold atoms are an ideal system to observe for the first 



time long photon storage in a system of N atoms in free 
space. Subradiance for two ions has been observed in 
the past |19j and a reduced decay rate into one radia- 
tion mode has been achieved for N atoms [50]. How- 
ever, it has not yet been possible to control and sup- 
press the decay into all vacuum modes for N atoms in 
free space extending thus the lifetime of the excitation to 
many times the natural lifetime of a single atom. Start- 
ing from the Ansatz used in previous work by several 
authors [TTJ |2"TH2~4] , we show that the exponential kernel 
of the dipolc-dipole coupling yields an important fraction 
of atoms to be coupled into subradiant modes. Following 
[2 5) , we study different inhomogeneous broadening mech- 
anisms, which allow to go beyond the fundamental Fano 
coupling by controlled coupling between super- and sub- 
radiant states using Doppler broadening and inhomoge- 
neous light shifts. In this letter, we consider three differ- 
ent experimental parameters to control subradiance: the 
optical thickness of the cloud, the driving laser intensity, 
and the temperature of the cloud. Moreover, we show 
how to distinguish subradiance from multiple scattering 
of light by tuning the driving laser frequency. 



We consider a Gaussian cloud, with root mean square 
size cr, of N two- level atoms (positions r iy transition 
wavelength A = 2ir/k, excited state lifetime 1/T), ex- 
cited by an incident laser (Rabi frequency Oqj detun- 
ing Ao, wavevector ko). We define the optical thickness 
6(A ) = b / [l + 4(A /T) 2 ] with b = 3A/(fccr) 2 its res- 
onant value. Restricting the atomic Hilbert space to the 
subspace spanned by the ground state of the atoms \G) = 
\g ■ ■ ■ g) and the single excited states \i) = \g ■ ■ ■ e 4 • • • g) 
and tracing over the photon degrees of freedom, one ob- 
tain an effective Hamiltonian describing the time evolu- 
tion of the atomic wavefunction \tjj) = a\G) + X^AN)- 
The effective Hamiltonian using standard approxima- 
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tions [2H can then be written as: 

HcS = [e lAot -' lko ^Sl +e- i * ot+ik °- T 'S , + ] 

i 

i i j^i 

where the first term describes the coupling to the laser 
field, the second accounts for the finite lifetime of the 
excited states, the third one describes the dipole-dipole 
interactions, with Vij = ^f^~^~p^ ; an d <S±) Si are 
the usual pseudo-spin operators for the kets \gi) and 
|ej). The effective Hamiltonian ([!]) describes dipole- 
dipole couplings in the scalar light approximation, where 
near field and polarization effects are neglected since we 
are considering dilute clouds, N(X/a) 3 <C 1. At low 
intensity (where the single excitation approximation is 
valid), a ~ 1 and the previous model describes a system 
of N classical dipoles driven by an incident electric field 
as expected by linear optics [27]- We use these classi- 
cal equations to study single photon subradiance. The 
driven steady state solution ~ \G) + e\TD) bears 
the phenomenon of single photon superradiance with a 
'timed Dicke' state \TD) = ^Ei eik ° r 'l«> and an am- 
plitude e ~ VNn /(2Ao + i(l + &o/12)r) 24J. This was 
e.g. exploited to explain the measured cooperative radi- 
ation pressure force in dilute clouds [23] . 

In order to study enhanced storage time based on Dicke 
subradiance, we will investigate the coupling between 
'short-lived' superradiant states, such as \TD), and the 
other states \<p) of the single excitation Hilbert subspace. 
Let us first consider the minimal coupling, based on ra- 
diative dipole-dipole coupling. As the effective Hamilto- 
nian Eq. is non-Hermitian, its eigenstates are non 
orthogonal ((^superlv'sub) 7^ 0) and subradiant states 
thus have common features with auto-ionizing states or 
Fano resonances [28]. This Fano-type coupling between 
'long-lived' subradiant states and 'short-lived' superradi- 
ant states leads to additional decay channels in addition 
to direct decay of the subradiant states to the ground 
state. This situation is reminiscent of the Hanle effect 
[29j where a competition of direct decay and transverse 
coupling can lead to surprisingly narrow resonances. For 
instance, the steady state solution of optical Bloch equa- 
tions for a system consisting of the ground state \G), one 
superradiant |y S uper) and a single subradiant state |iyJ S ub) 
reveals that in the absence of Fano coupling the steady 
state solution of the subradiant state is zero. Neglecting 
the direct decay of the subradiant state to the ground 
state on the other hand, one finds that for resonant ex- 
citation and small coupling between the excited states, 
after a long transient time all atoms are pumped into the 
subradiant state |y su b)- The optical Bloch equations of 
such a simplified three-level model show that, by increas- 
ing in a controlled manner the coupling between super- 
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FIG. 1: (color online). Time evolution of the normalized ex- 
cited state population (black solid curve) after switching off 
the laser for N = 2000 atoms, ka = 10, A = 10 T (the 
laser was on before during 50 T -1 to let the system reach the 
steady state). At first, the population decreases faster than 
the single atom decay (black dashed line) and then slower. 
We respectively identify these phenomena as super- and sub- 
radiance. The inset shows the emission diagrams of the su- 
perradiant 'timed Dicke' state \TD) (blue) and subradiant 
modes (red) averaged over 8 realizations (rescaled to allow 
convenient comparison) . 



and subradiant states, it should be possible to efficiently 
store populations into 'long-lived' subradiant modes. 

In Fig. [I] we show the normalized excited state popu- 
lation cx \Pi\ 2 as computed from numerical solution of 
the effective Hamiltonian Eq. ([I]) in the linear regime. As 
precise initial conditions play a crucial role in the subse- 
quent fast and slow decay [27] , we start initially with all 
atoms in the ground state \G) and keep the coherent laser 
drive for 50 T -1 before switching it off, realizing thus ex- 
perimentally accessible conditions. The fast initial decay 
of the superradiant state r supcr = (1 + b /12) T is clearly 
seen. Moreover, after this initial fast decay, subradiance 
manifests itself in a slowly decaying excited population 
with a rate well below the single atom decay rate. At 
first, the subradiant decay is not purely exponential since 
several modes decay simultaneously. For longer times, it 
then ends up with a pure exponential decay (referred as 
subradiant decay in the following) when only one 'long- 
lived' mode dominates. The emission diagram of the su- 
perradiant 'timed Dicke' state \TD) is clearly forward 
directed (see Fig. [I]), a phenomenon reminiscent of Mie 
scattering. On the other hand, subradiant modes show 
isotropic diagrams. They do not possess the symmetry of 
the laser excitation since they are not directly coupled to 
it, a feature which can be exploited in the experimental 
detection of subradiance. 

In Fig. [2j the subradiant decay rates (due to Fano 
coupling) are plotted as a function of the inverse on- 
resonance optical thickness of the system. The difference 
to multiple scattering of light which can also yield long 
photon trapping times, can be understood by looking at 
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FIG. 2: (color online). Subradiant decay rate as a function of 
l/6o (we kept N = 400 constant, Ao = 10 F and varied ka). 
At large optical thickness it scales as oc F/feo (black dashed 
line) and saturates to V for dilute clouds. Fhe insert shows 
the subradiant decay rate as a function of detuning for b = 3 
(N — 400, ka — 20). Fhe shaded area corresponds to the 
multiple scattering region b(Ao) > 1. 



the decay rates for different excitation frequencies. As 
shown in the inset of Fig. [2] close to resonance we ob- 
serve reduced decay rates, which we associate to the large 
optical thickness for resonant photons [5]. For larger de- 
tunings however, the decay rate becomes independent of 
the excitation frequency, consistent with the subradiant 
nature of states weakly excited off-resonance. 

The fundamental Fano coupling between sub- and su- 
perradiant states in the Dicke basis can be understood 
from the diagonal terms in the bare bases \i). Indeed, 
the local field at i-th atom is the sum of the external 
field Eq and the field scattered by all the other dipoles j 
at the location of the atom i: 



E tot (r t ) = Eo^+J^E^i) - £o(r;)(l 



■ £i.e 



')■ (2) 



We describe the scattered field by a small local speckle 
field with a random amplitude scaling as Si cx ^fe(Ao) <C 
1 and a random phase ifi . The amplitudes Pi depend on 
the local field and differ from the driven timed Dicke am- 
plitudes. The dipole i is thus driven with a random field 
corresponding to an inhomogeneous broadening mecha- 
nism |25) . where we now have inhomogeneity in ampli- 
tude and phase. An estimation of the expected lifetimes 
of the 'long-lived' subradiant modes is also an important 
issue. We have checked numerically they scale as 
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for 6o > 1 (see Fig. [2]). This scaling is close to what can 
be obtained assuming that in the limit of large detun- 
ing Ao, where multiple scattering inside the sample can 
be neglected, the escape rate of the excitation from the 
sample is well approximated by the inverse of the 'long- 
lived' mode lifetime, which scale as 1/bo (30] . This can 



FIG. 3: (color online). Super- (blue curve) and subradiant 
(red curve) decay rates as a function of Doppler broadening 
for N = 200 atoms, ka — 5 and a laser detuning of Ao = 10 F. 
Initially the laser is driving the atoms for 50 F -1 and the 
subradiant decay rate is evaluated 50F -1 after the laser is 
switched off. 



also be compared to scaling laws obtained from quantum 
chaotic scattering theory [3T] : for large bo (corresponding 
to a large number of atoms N compared to the number of 
outgoing modes M ~ {ka) 2 ) the minimum width of the 
resonance is expected to scale as r sub ~ M/N ~ 1/bo. 

We now turn to the possibility of controlled coupling 
between super- and subradiant states, opened by cold 
atoms. In atomic physics many different inhomogeneous 
broadening mechanisms are known. We will focus on 
two such mechanisms, well adapted to become a control 
knob to steer excitations into the subradiant state. Fet 
us thus consider the impact of residual motion of the 
atoms. The effective Hamiltonian can be extended to 
include Doppler shifts and time dependent positions of 
the atoms. The coupled equations for the dipole ampli- 
tudes Pi read in the linear regime 
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where Vij(t) 



fclri-r-f-Kvi-Vj)*! 

Solving these equations for increasing temperature, we 
notice that the fast superradiant and the slow subradi- 
ant decay rates draw closer to the single atom decay rate 
r (see Fig. [3]). This result shows that the fragile subradi- 
ant modes are quickly destroyed even by moderate atomic 
motion. Cold atoms seem a well adapted system allow- 
ing to tune the atomic motion from being negligible to 
becoming dominant. We also checked that the dominant 
term in the reduced super and subradiance stems from 
the position dependent dipole-dipole coupling Vij (t) term 
rather than from the random detuning term in Eq. Q. 
This dependence on the atomic motion explains why sub- 
radiance of N atoms has not been observed in hot atomic 
vapors, despite the efforts in this field in the 1970s [5]. 
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Exploiting the possibility of controlled coupling via in- 
homogeneous broadening even further, we studied the 
role of larger laser intensities on super- and subradiance. 
As the effective Hamiltonian approach is only valid to 
first order in the Rabi field coupling, it can only treat the 
linear optics regime. We therefore used a master equa- 
tion approach, where the evolution of the atomic den- 
sity operator p in the electric-dipole, rotating-wave, and 
Born-Markov approximations is given in the interaction 
picture by [521 155] 

p = ^ (h cSP p hI s ) ytjsLpSX (5) 

* 3 

where in the scalar light limit 7^ = r s '"|^. _ r rp . Pro- 
jecting Eq. ([5]) on the different Fock states, we obtain 
a set of coupled equations for the density matrix ele- 
ments p G \ G = (G\p\G), Pi\ G ee (i\p\G) and p t y ee (i\p\j). 
This approach is still restricted to the same Hilbert sub- 
space with at most one excitation and is thus limited to 
moderate laser intensities or to situations where multi- 
ple excitations are suppressed as for instance in the case 
of Rydberg blockade. However, it does not require the 
ground state population to remain unaffected and can 
take into account the light shifts of the states. We have 
checked analytically and numerically that at first order in 
Rabi frequency of the laser coupling, the master equation 
is equivalent to the equation used in the effective Hamil- 
tonian approach Eq. ([I]) in the linear regime a ~ 1 (with 
correspondence Pg\g 1; Pi\G Pi)- Here we exploit 
the small fluctuations induced by the random local field 
driving the individual atoms i. As the local field has a 
random speckle structure (see Eq. Q), the random light 
shifts and phases can be understood as an inhomoge- 
neous broadening mechanism, depending on the interfer- 
ence term between the incident field £t)( r i) and the scat- 
tered field Eie lVi Eo(ri). The advantage of the light shift 
coupling and dephasing is the flexibility it offers as the 
laser intensity can be easily and quickly controlled. Fig. 
[4] shows the super- and subradiant decay rates as a func- 
tion of the Rabi frequency Slo and the subradiant fraction 
(i.e. the remaining excited state population) after 
switching off the laser (the laser was on during 50 T -1 ). 
When the Rabi frequency was varied we checked that the 
excited state population remains small (J2iPi\i < 0.15) 
to ensure consistency with the model. We observed in 
Fig. [4] (b) up to a three time increase of the subradi- 
ant fraction when the intensity is raised compared to the 
£lo — > value determined by Fano coupling. Note that 
the subradiant fraction would be the same for any inten- 
sity with the effective Hamiltonian approach (due to lin- 
earity). Fig. [4]illustrates how changing laser intensity al- 
lows controlling the coupling strength between super- and 
subradiant modes and subsequently the subradiant pop- 
ulation. As the laser is switched off, the local field, based 
on the interference between the incident and scattered 
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FIG. 4: (color online), (a) Super- (blue curve) and subradiant 
(red curve) decay rates as a function of laser intensity for 
N — 200 atoms, ka — 5 (60 = 24) and a laser detuning of 
Ao = 10 r. (b) Subradiant fraction for the same parameters. 
The laser is switched on during 50 T -1 and the subradiant 
decay rate and subradiant fraction (i.e. the remaining excited 
state population) are computed 50r _1 after switching off the 
laser. 



field, is quickly reduced by a large amount. The inho- 
mogeneous coupling is thus significantly decreased clos- 
ing the 'door' between the subradiant and superradiant 
states. One can see this effect in Fig. [4] (a), where a de- 
crease in the superradiant decay rate is observed because 
there is still some inhomogeneous broadening source just 
after switching off the laser. However, for longer time the 
local field is much less intense and almost no variation of 
the subradiant decay rate is seen - light induced inhomo- 
geneous broadening source is no longer present. In that 
case the subradiant decay rate is just the same as the 
one given by Fano coupling. In the same way, for low in- 
tensities, the inhomogeneous broadening induced by the 
laser is dominated by the Fano coupling. The subradi- 
ant decay and subradiant fraction then remain almost 
unaffected, as illustrated in Fig. [4] for J7o < 0.1 T. 

In conclusion, we have shown that inhomogeneous 
broadening schemes allow to understand and control stor- 
age of an optical excitation into 'long-lived' subradiant 
modes. We have proposed to use the cloud optical thick- 
ness, the driving laser intensity, or the cloud temperature 
as possible experimental control parameters for subra- 
diance, but further parameters, such as a far detuned 
speckle field or magnetic field would yield similar results. 
This opens the door for the first observation of Dicke 
subradiance of photons in a cold cloud of N atoms in 
free space. Inhomogeneous broadening schemes will also 
be of interest to studies of Anderson localization of light 
in resonant two-level systems [ID] . Mapping the inho- 
mogeneous coupling schemes described in this letter to 
a lambda scheme as used in quantum information sci- 
ence [51] will help to understand limitations of storage 
of qubits in atomic vapors in modes less exposed to fast 
decay in a decoherence free subspace [35]. Controlled 
transfer to other modes can also be achieved with non- 
random coupling, as e.g. used in slow light experiments 
[36] . Thus we expect that novel schemes to engineer more 
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robust and faster storage and exploit the larger Hilbcrt 
space can be addressed using two-level systems as toy 
models. 
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